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Abstract 

We examine the long-term time-dependence of Gaussian wave packets in a circular infinite well 
(billiard) system and find that there are approximate revivals. For the special case of purely m = 
states (central wave packets with no momentum) the revival time is Tr'^~'^^ = Sfj.R'^/hTT, where 
fj, is the mass of the particle, and the revivals are almost exact. For all other wave packets, we 
find that Trev^^'^ = (tt"^ /2)Tr^~^^ ~ 5Trev~'^^ and the nature of the revivals becomes increasingly 
approximate as the average angular momentum or number of m 7^ states is/are increased. The 
dependence of the revival structure on the initial position, energy, and angular momentum of the 
wave packet and the connection to the energy spectrum is discussed in detail. The results are also 
compared to two other highly symmetrical 2D infinite well geometries with exact revivals, namely 
the square and equilateral triangle billiards. We also show explicitly how the classical periodicity 
for closed orbits in a circular billiard arises from the energy eigenvalue spectrum, using a WKB 
analysis. 
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I. INTRODUCTION 



The connections between the quantized energy eigenvalue spectrum of a bound state and 
the classical motions of the corresponding classical system have become increasingly interest- 
ing and important as the ability to experimentally probe the quantum-classical interface has 
dramatically improved. Methods such as periodic orbit theory , for example, provide 

direct connections between the energy spectrum and the closed orbits of the classical system. 
The study of the time-dependence of wave packet solutions also depends critically on the 
energy spectrum, most especially as related to the existence of revivals (and superrevivals), 
wherein initially localized states which have a short-term, quasi-classical time evolution, can 
spread significantly over several orbits, only to reform later in the form of a quantum revival 
in which the spreading reverses itself, the wave packet relocalizes, and the semi-classical 
periodicity is once again evident. Such revival phenomena have been observed in a wide 
variety of physical systems, especially in Rydberg atoms P], and calculations exist for many 
other systems 

The archetype of a one-dimensional model system for quantum revivals is the infinite well 
(where such revivals are exact) and a number of analyses p - have provided insight into 
both the short-term and long-term behavior of wave packets. Just as with many systems 
which depend on a single quantum number, one typically expands the energy eigenvalues 
(assuming integral values) about the central value used in the construction of a wave packet 
via 

E{n) ^ E{no) + E'{n^){n - no) + ^^"(no)(n - no)^ + \E"\n^){n - no)^ + ■ ■ ■ (1) 

/ o 

in terms of which the classical period, revival, and superrevival times are given respectively 
by 

^''-W¥o)\ PpfP |E-(no)|/6 

Systems with two quantum numbers with energies labeled by E{ni,n2), offer 

richer possibilities for wave packet revivals and typically the long-term revival structure 

depends on three possible times, given by 

rim) ^ 27rn ^ 27rn 

{l/2)d^E{ni,n2)/dnl {l/2)d^E{ni,n2)/ dnl ^ ' 

and 

d^E{nun2)/dnidn2 ^ ' 



and the revival structure depends on the interplay between these three times. 

The two-dimensional generalization of the infinite well, or two-dimensional infinite square 
well (billiard) of size L x L, provides the simplest example of such a two quantum number 
system, and in this case the revival times are identical, namely 

rev rev \^/ 

with no cross-term present. For rectangular infinite wells with incommensurate sides [L^ x 
Ly, Lx/Ly 7^ p/g), the structure of the revival times may be more complex 12|. 



Another trivially related 2D billiard system for which exact quantum revivals are guar- 
anteed is an isosceles triangle infinite well with 45° angles. The solutions for that system 
can be obtained from the 2D square well (cut along the diagonal) by taking appropriate 
linear combinations of the form [u{nx){x)u{ny){y) — U{n^){y)u{^ny){x)]/ \/2 which then sat- 
isfy the boundary conditions along the long side (hypotenuse) of the triangle, so long as 
nx ^ Uy. The energy spectrum therefore consists of one set of the 2D square well energies, 
E{nx,ny) = h^7!''^{nl + n^)/2mL^, but restricted to n^. < Uy. The purely quadratic energy 
dependence guarantees that all wave packets will have the standard revival time given by 
Eqn. ©. 



A less obvious case of exact quantum revivals in a novel 2D bif 



iard shape exists for the 



14| (with side of length L) 



equilateral triangle. The energy spectrum for this system jl^l - 
is given by 

E{p,q)=(^^^ (^-^^^ (P^ + (1^ - (IP) withl<g<p/2 (6) 

with all states being degenerate (due to symmetries about obvious axes), except for those 
with p = 2q. Because of the exact quadratic dependence (with trivially related coefficients) 
on the two quantum numbers, there are exact revivals with a common revival time given by 

rev rev rev 4/),7r 

Given the fact that these two simple geometrical cases of 2D infinite well (billiard) sys- 
tems exhibit exact quantum revivals, it is of interest to study to what extent, if any , the 
circular infinite well will exhibit revival structures. Just as square/rectangular and 
spherical/circular billiard geometries were among the first considered using periodic or- 
bit theory, a comparison of the revival structure in these two distinct systems is appropriate. 



While qualitatively different from the equilateral triangle or square wells, one can perhaps 
suggestively describe these geometries as regular A^-sided polygons with = 3 and = 4 
respectively, so that the circular case we will consider here corresponds to = oo. 

The ability to fabricate such billiard systems (or analogs thereof) and experimentally 
probe the energy spectrum or time-dependence makes such a study of more than academic 



interest. For example, the energy- level structure and statistics of microwave cavities [17 1 
provide an analog system to probe periodic orbit theory and statistical measures of chaotic 
behavior in arbitrary shaped 2D billiard geometries. Measurements of conductance fluc- 
tuations in baUistic microstructures [3] have been tentatively used to identify frequency 
features in the power spectrum with specific closed orbits in a circular (and stadium) bil- 
liard. More recently, the realization of atom-optics billiards 19], with ultra-cold atoms in 
arbitrary shaped 2D boundaries confined by optical dipole potentials, has allowed the study 
of various chaotic and integrable shapes such as the stadium, ellipse and circle, again for 
short-term, semi-classical propagation. Motivated by these studies, we will focus on the 
long-term revival structure of quantum wave packets in a circular billiard geometry. 

Since the revival structure of any system depends on the quantum number dependence 
of the quantized energy eigenvalues, in Sec. II we discuss, in some detail, the dependence 
of the two-dimensional circular infinite well eigenvalues on the radial and angular momen- 
tum quantum numbers, rir and m. We then, in Sec. Ill, present numerical results for the 
autocorrelation function for several classes of initially Gaussian wave packets in the circular 
well. 

Based on both the energy eigenvalue analysis and our detailed wave packet simulations, for 
purely m = wave packets (implying radially symmetric states, with no average momentum) 
we find almost exact revivals with Tr™~°'' = SfiR'^/hn where /i denotes the mass of the point 
particle in the circular billiard of radius R. For wave packets including |m| 7^ components, 
meaning any with non-vanishing average momentum or not initially localized at the center 
of the well, there are only approximate revivals, becoming increasingly so as the angular 
momentum is increased. The revival times in all these cases is Tr^^^"* = {tt^ /2)Tr'^~^^ 
5TJI^5~°^ due to the seemingly accidental fact that IO/tt^ = 1.013. 

Finally, in one Appendix (A) we point out a similarity between a pattern of special 
'accidentally' revival times in the 2D circular, square, and equilateral triangle wells, while 
in a second Appendix (B), we show explicitly how the classical periodicity for closed orbits 



4 



in the circular well arises directly from the quantum mechanical energy spectrum, using a 
WKB analysis. 



II. ENERGY SPECTRUM FOR THE CIRCULAR INFINITE WELL 

The problem of a point particle (with mass denoted by fi, to avoid confusion with the 
angular momentum quantum number, m) confined to a circular infinite well of radius R is 
defined by the potential 

{0 for r < 
(8) 
oo for r > R 

The (unnormalized) solutions of the corresponding 2D Schrodinger equation are given by 

i:{r,e) = J\m\{kr)e'"'' (9) 

where the quantized angular momentum values are given by = mh for m = 0, ±1, ±2, ... 
and the J\m\{kr) are the (regular) Bessel functions of order \m\. 

The wavenumber, k, is related to the energy via k = ^J^^lEJW and the energy eigenvalues 
are quantized by application of the boundary conditions at the infinite wall at r = i?, namely 
J\m\{z = kR) = 0. The quantized energies are then given by 

Em = 2;.i?2 (10) 

where Z(^m,nr) denotes the zeros of the Bessel function of order \m\, and counts the number 
of radial nodes. The energy spectrum is doubly degenerate for \m\ ^ corresponding to 
the equivalence of clockwise and counter-clockwise (m > and m < 0) motion. Because 
the quantum number dependence of the energy eigenvalues is the determining factor in 
the structure of wave packet revivals, we wish to examine the m, Ur dependence of the 
E(m,nr) a: [^{m,n,)]^ m detail. 

As a first approximation, one can look at the large z behavior of the Bessel function 
solutions [2^ for fixed values of |m|, namely 

■'hW — i/J<=o«(^-^-^)+--- (11) 

With this approximation, the zeros would be given by 

\m\ vr TT / |m| 3\ 

Zim,nr) 2 4~(2'^'- + l)2 °^ Z(m,nr) ^ \nr + — + -j ir = ZQ{m,nr) (12) 



where we define the function zo{m,nr) for future reference. If this result were exact, the 
quantized energies would depend on two quantum numbers in at most a quadratic manner 
and there would be exact wave packet revivals, just as for the 2D square or equilateral 
triangle billiards. 

To see to what extent this approximation is, in fact, valid for various values of m, n^, 
we evaluate (numerically) a large number of the lowest-lying 'exact' Bessel function zeroes, 
Z{m,nr)- We then plot, in Fig. 1, the combination 2z(m„^)/7r — 3/2 versus m which would 
yield horizontal lines corresponding to constant, integral values of 2nr + \'m\ if the result 
in Eqn. (fT^ were exact. For three such states, we also illustrate the corresponding radial 
probability density and note that for increasing values of |m|, there are fewer radial nodes 
so that the appropriate values oi z = kR used in the boundary condition are smaller and 
the approximation used in Eqn. (fT^ becomes worse. 

We see that the approximation is only a good one for |m| with obvious quadratic 
corrections. Instead of attempting to evaluate Bessel function zeros to higher precision us- 
ing more elaborate 'handbook' expansions of J\m\{z) for large z, we note that the result of 
Eqn. fll2|l . and the important and necessary corrections to it, can be obtained in a straight- 
forward and accessible way by use of the WKB approximation. 

If we first quantize the angular variable to find that angular momentum is given by 
Lz = mh, we can note that in the radial direction the particle moves freely up to the 
infinite wall at r = i?, but is subject to an effective centrifugal potential given by Ve//(r) = 
L^/2yur^ = (m^)^/2/ir^. The classical particle cannot penetrate this centrifugal barrier and 
has an associated inner radius (distance of closest approach) given by 

, m^/j,^ \m\h , . 

VeffiRmin) = = ^ -^mm = —J^=^ (13) 

We can also write this in the useful form 

^ \m\R , ^ h^z"^ 

Rmm = ^ Where - ^ (14) 

more directly in terms of the desired dimensionless variable which is equivalent to the 
energy eigenvalue. In Fig. 1, we have also indicated the value of Rmin/R = \m\/z for several 
of the radial probability densities, to compare to the fully quantum mechanical results. 
The WKB quantization condition on the radial variable, r, is then given by 



kr{r)dr={nr + CL + CR)7T where K{r) = \l I - ^ (15) 



rlR(r)P vs. r 



N 




r=0 R . /R r=R 

min 

FIG. 1: Plot of the lowest-lying Z(m,,nr) zeroes of the Bessel function, determined (numerically) by 
J\m\{z{m,nr)) = versus TTi, scalcd as 2z(^ „^')/7r — 3/2. Zeroes corresponding to fixed values of the 
combination 2nr + \m\ are shown connected, not 'fit' to data. Shown on the right are plots of the 
radial probability density, P{r) = r|i?(r)p versus r for the three states shown with a 'circle' for 
which a value of 2nr + \m\ = 20 is fixed. The value of Rmin/R = \m\/z suggested by Eqn. (|14|) is 
also indicated by a vertical dashed line. 



The matching coefficients j2l| are given hy Cl = 1/4 and Cfi = 1/2 which are appropriate for 
'hnear' boundaries (at the inner centrifugal barrier) and 'hard' or 'infinite wall' boundaries 
(such as at r = i?), respectively. The WKB energy quantization condition for the quantized 
energies, in terms of n,. explicitly and \m\ implicitly, via both the E and Rmin terms, can 
then be written in the form 



{rir + 3/4)7r 



(16) 
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The integral on the left can be evaluated in the form 



™" dr = Jr^ - - sec-i -- ] (17) 



Rri 



Y mm -"//tin--- 1 p 



R 



\J\ — — X sec ""^(l/a;) 



where we define x = Rmin/R = \'m\/z. This result can be expanded for small values of x 
(i.e., Rmin/R « 1 or \m\/z « 1) to obtain 



TTX x" Sx*^ 



Vl -x2 -xsec^Vl/x) = 1 + — + — + — + + ■•• (18) 

^ ' ' 2 2 24 80 896 ^ ' 

The WKB quantization condition in Eqn. |TB|) can then be written, in terms of z, in the 
form 

If we keep only the first two terms on the left-hand side, we find that 

z (nr + \m\/2 + 'i/A)'K = zoijn.nr) (20) 

which is the lowest-order result obtained directly from the limiting form of the wavefunction. 

To improve on this result, we simply keep successively higher order terms, solving itera- 
tively at each level of approximation using a lower-order result for z, and we find the much 
improved approximation 

_ 7 83 6949m8 

Z{m,nr) - Zo[m, Ur) - 2^^^^^ ~ 24 [2o(m, n^)]^ ~ 240 [zo{m,nr)f " 13AAO[zo{m,nr)Y 

(21) 

which we have confirmed numerically is an increasingly good approximation, especially for 
Ur >> 1. For the study of wave packet revivals, we only require the energy eigenvalue 
dependence on m, to second order, but higher order terms such as those above might be 
useful for super-revivals and even longer-term time-dependence studies (or more detailed 
analytic periodic orbit theory studies of the circular well.) 

For the special case of m = 0, we find no improvement using this WKB technique, but 
motivated by the form of the expansion in Eqn. (j2^, we fit the first 50 lowest-lying m = 
zeros to a similar form and find the result 

^,„,„,, = ^„(0, ./.) + - 24i,„,J_„^)p + ■ • • (22) 
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We cannot unambiguously fit to any higher-order terms, as much of the non-hnear spacing 
information is contained in the lowest few zeros. 

Using Eqns. (j2ip and (j22j) . we can evaluate the energy eigenvalues to quadratic order in 
rir, m in order to probe the revival structure of wave packets. For the special case of m = 0, 
we find that 



{0,nr) 



2/ii?2 2/ii?2 

while for the more general case with m 7^ 0, we find 



3\' 1 
4 / 47r2 



(23) 



(m,nr) 



2/ii?2 2/ii?2 



m 

rir H h 

" 2 4 



3 \ m 



(24) 



The fact that these energies depend on non-integral values of the effective quantum 
numbers is reminiscent of the case of Rydberg wave packets in alkali-metal atoms due to 



quantum defects 



22| and methods similar to those used there might prove useful. In what 



follows, however, we simply examine the time-dependence of typical m = and m 7^ wave 
packets directly. 

III. GAUSSIAN WAVE PACKETS AND REVIVALS 

Any wave packet in the circular billiard can be expanded in the normalized eigenstates 
of the form 



where 



R 



[N{ra,nr)\ / r [j\m\{kr)\ dr = 1 
Jo 



with expansion coefficients given by 



Hm,nr) = {'ip{r,0;t = 0)\lp(m,nr)) 



which satisfy 



+00 00 



^{m,nr) I 



1 



(25) 
(26) 

(27) 
(28) 



m=— 00 nr=0 

The expectation value of the energy in this potential well is given by 



m=—oonr=0 ^ ' 



(29) 



and expectation values of powers of angular momentum are also easily evaluated to give 



-oo oo 



(30) 



m=—oo nr=0 

The subsequent time- dependence of the wave packet is then given by 



-oo oo 



m=—oo 71^=0 

and the standard autocorrelation function [23] is given by 



(31) 



hoo oo 



m=—OD nr=0 

For definiteness, we will use a standard Gaussian wave packet of the form 



(32) 



where 



'4!o{x;xo,pox,b) 



^ipOx{x-xo)/h^-(x-xo)'^ /2b'^ 



(33) 



(34) 



with a similar expression for ipo{y; yo, poy, h) . The initial expectation values for the x variables 
are given by 



and 



(x)o = Xq 



{Px)o= POx 



(x^)o = xl + 
(pDo = pIx + 



2 



Axn 



V2 



n 



(35) 



(36) 



262 ' 

with similar results for y. So long as the initial location, {xo,yo), is well away from the 
edges of the potential well, such a Gaussian form can be easily and reproducibly expanded 
in terms of eigenstates. The expectation value of total energy is 



1 

2m 



21 



{POxf + {PQyf + 



62 



(37) 



In this central potential, angular momentum is conserved and we also have the specific 
results for this Gaussian form 



and 



(L^) = {xpy - ypx) = {x){py) - {y){Px) = xoPoy - yoPox 



62 ^2 

(L2) = (Xopoy - yoPoxY + ^ [{Poxf + iPOyf] + ^ [(^o)' + (l/o)'] 



(38) 



(39) 



10 



so that 



{Am)h = AL, 



62 



[{Poxy + {p 



(40) 



As a check on the numerical evaluation of the expansion coefficients, it is useful to be able 
to compare the general results for (E) and {Lz ) in Eqns. fl29j) and with the specific 
results for the Gaussian in Eqns. (jSZj)? (ISHl), and (jH^jl . 

We begin by focusing on the special case of zero-momentum wave packets centered at 
the origin, namely with vanishing values of {pox,Poy) and (a;o,?/o) in which case the initial 
wave packet is radially symmetric and therefore has an expansion in pure m = angular 
momentum states. (This is consistent with the result in Eqn. pn|l which has AL = Am = 
for this state). 

For such states, where only the m = eigenstates contribute, we can write the energy 
eigenvalues from Eqn. (|^^ in the form 

2 1 



Eirir 



32/ii?2 



n,. + 



+ 



An' 



O 



[rir + 3/4y 



Iriri^rir + 3) + 

N 1 
Knr) + o + 



9 + 



(41) 



8 27r2 



where /(n,.) = nr{2nr + 3) is an integer (neither even nor odd in general). The last term 
in the square brackets is independent of rir and will make the same, constant, overall phase 
contribution to the autocorrelation function, so we focus on the /(n^) term. Since this integer 
has no special evenness/oddness properties, its contribution to the phase of each |a(„^„^)|2 
term in Eqn. (jH^ will be identically unity at a revival time given by 



rp{m=0) 
J rev 

h 



2n 



or 



rp{m=0) _ ^ 



2^i?2 



4Tn 



(42) 



V4/ii?2 

Thus, at integral multiples of 4To, we expect nearly perfect revivals because of the almost 
regularly spaced structure of the m = Bessel function zeros. At these recurrences, we can 
also predict the overall phase corresponding to the last term in Eqn. (PT|) . namely 



-ift2^2/4M-R^(4To)(9/8+l/27r2) 



-27ri(9/8+l/27r2) 



g-27rig-27ri(l/8+l/27r2) ^ ^-ittF ^^g-j 



where F = 1/4 + l/vr^ ^ 0.351. 
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To investigate these predictions numerically, we have used a Gaussian of the form in 
Eqn. with the specific values 

1 



2m = h = R = l 



and 



IOV2 



so that Axo = A?/o 



0.05 



(44) 



Using the normalized eigenstates, we numerically evaluate the overlap integrals to obtain the 
ci'im,nr)y using enough states to ensure that the appropriate conditions, such as Eqns. ((23), 
in^ . and (jHUIl . are all satisfied to better than 10~^ accuracy. 

Using the expansion coefficients for this state, we plot the modulus squared of the auto- 
correlation function, |y4(t)p, in the bottom plots of both Figs. 2 and 3, with time 'measured' 
in units of Tq. The almost exact revival structure at integral multiples of 4Tq is evident. 
As a further check, we can evaluate the phase of A{t) at each revival and find that to an 
excellent approximation it is given by —nFir as in Eqn. (j43j) . If one decreases (increases) 
the value of b, so that the initial wave packet is narrower (wider), the energy eigenvalues 
required to construct the packet are then larger (smaller) (from Eqn. (j37j) ) and are therefore 
generally more (less) evenly spaced (from Eqn. ((211)) and we indeed confirm this with our 
numerical simulations; the eventual, long-term decrease |y4(nTo)| with increasing n is faster 
(slower) for smaller (larger) values of b. 

We next move away from the special case of the zero-momentum, central wave packet by 
considering individually the case of Xq 7^ and pQy 7^ (but not both). In each case, the 
average angular momentum of the state is still vanishing (from Eqn. (jHHj) ). but m 7^ values 
of the expansion coefficients are now required. We must now use the more general case for 
the energies, which to second order in m 7^ 0,nr, are given by Eqn. 

2 



TT 



(m,nr) 



32/ii?2 
32/ii?2 



\m\ 



m 



[I6nl + 2Anr + 16\m\nr) + 4|m|(|m| + 3) - 



16777.2 



(45) 



^rir) + 8l{\m\) 



16m2 



+ 9 



where 



l[nr) = nr{2nr + 3 + 2\m\) 



and 



/(|m|) = \m\{\m\ +3)/2 



(46) 



are both integers, again, with no special even or oddness properties. We can then write 
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P0y=0 



t/T„ 



x„/R 
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FIG. 2: Plot of the autocorrelation function, versus t, in units of Tq = 2/ii?^/^7r, The 

numerical values of Eqn. ()44() are used, along with yo = and pox = Poy = 0. The results for 
versus t, as one varies the xq of the initial wave packet away from the center of the circular 
billiard, are shown on horizontal axes. 



these energies in the form 



(m,nr) 



27ih 



i{nr) + li\Tn\) - ^ + ? 



(47) 



At integral multiples of the m = revival time, = N{4To), the first two terms give 
g-Af(27rj) _ ^ phases to each (m, rir) term in the autocorrelation function, while the last term 
gives an overall, (m^n,,)-independent phase, just as in the m = case. The other term, 
however, gives a contribution 



-(2ni){m^N){2/ 



(48) 



which depends on m explicitly and which therefore eliminates the revivals, increasingly so, 
as the wave packet is dominated by m 7^ terms. However, because of a seeming numerical 
accident, at integral multiples of 5Tr^~°'' = 20To, we recover approximate revivals due to 
the fact that 5 x (2/7r^) = 1.013. We thus find approximate revivals for the more general 
m ^ case given by Tje"^°^ = (7rV2)T4'S=°) ^ 5X4^=^°)^. 

This effect is illustrated in more detail in Figs. 2 and 3 where we plot |^(t)P versus t 
as we move from the central, zero-momentum wave packet by first moving away from the 
origin (xq 7^ in Fig. 2) or having non-zero momentum values {poy 7^ in Fig. 3). In 
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FIG. 3: Same as Fig. 2, but with xq = yo = and pox = 0, as one increases poy. 

each case, as we increase the parameter (xq or Poy), we necessarily include more and more 
|m| 7^ eigenstates. For even a small mix of such states, the Tr'^~'^^ revival periods at 
most integral multiples of ATq disappear, while evidence for the more general T^eif^^ = 20To 
revivals remains evident. 

For the particular numerical values used in Eqn. the spread in angular momentum 
required from Eqn. pn|) is given by 



AL 



Am 



(49) 



(since h = 1) so that the x 7^ and poy 7^ values used in Figs. 2 and 3 actually correspond 
to the same set of AL for each horizontal case shown. 

We note that this 'lifting' of a seemingly 'accidental' degeneracy in the pattern of revival 
times is somewhat similar to the special case of a zero-momentum Gaussian wave packet 
in a 2D square or triangular billiard, initially placed at the center, cases which we briefly 
discuss in Appendix A. 

This pattern of revival times depending on two distinct quantum numbers is also some- 
what reminiscent of that encountered in a rectangular billiard with differing sides of length 
Lx,Ly where if the sides are incommensurate one would expect a less elaborate revival 
structure. Since the revival times typically scale as Trev oc L^, the appearance of a vr^ scale 
factor which can give rise to very close to an integer ratio lO/vr^ ^ 1 (to within 1.3%) 
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is appropriate; in this case, the relevant length scales for the radial quantum number and 
azimuthal quantum numbers are most likely multiples of R and 27iR respectively, so that 
relative factors of vr^ in the revival times can appear naturally. 

The presence of the Am 7^ revivals becomes increasingly less obvious as the average 
angular momentum is increased away from zero (with both Xq and poy now non- vanishing) , 
since the required energy eigenvalues are in a region of large \m\/z where the lowest-order 
approximation (from Eqn. ((201)) of evenly spaced z values becomes worse. We also note that 
even with (L) = 0, as we increase p^y, the spread in m values required also increases (as in 
Eqn. pn|l ). so that the overall number of states required to reproduce the initial Gaussian, 
and which have to 'beat' against each other appropriately, increases as well, making revivals 
more difficult to produce. The increasingly large number of states required to construct the 
Gaussian wave packets for larger values of (L) can also be seen during the collapsed phase 
when the 'average value' of |A(t)p, namely X^mn^ k(m,nr) 1"^; becomes increasingly small as 
the fixed probability (constrained via ^^^^ |fl(m,,nr)P = 1) is spread over more and more 
states. 

While we have focused on the long-term, revival structure of the autocorrelation function, 
the appearance of more short-time features in |A(t)|, corresponding to short-term, semi- 
classical closed orbits, is also apparent. Such trajectories are characterized [l^ by periodic 
orbits with path lengths and minimum radii (distances of closest approach) given by 

L(p, q) = 2pRsm ( — ) and Rmin = Rcos ( — ) (50) 
\p J \p J 

with integral values of {p, q) (with p > 2q) describing the number of 'hits' on the 
wall and the number of 'revolutions' for one complete orbit respectively. If, for exam- 
ple, we place Gaussian wave packets with = and poy > at locations given by 
(xo,?/o) = {Rmin{p,(l),0), "we find obvious peaks in the autocorrelation function at times 
given by Tci{p,q) = L{j),q) / {jpoy/m) corresponding to classical closed orbits. (This struc- 
ture is evident, of course, only when the expected classical periods are less than the wave 
packet spreading time. At; this can be estimated using the result for a free Gaussian as 
At = 2mAxl/h.) We discuss in Appendix B exactly how the classical closed orbit pe- 
riodicity is reproduced from the quantum mechanical energy spectrum, using the WKB 
approximation of Eqn. 

Variations on the problem of a circular infinite well can also be examined for their possible 
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revival structure. The 'half-circular' well, with an infinite wall added along a diameter, is 
exactly soluble with linear combinations of e*"*^ and e"*"*^ solutions being able to satisfy 
the new boundary condition for m 7^ 0, while the m = solutions are not allowed. The 
energy spectrum then consists of one copy of the -E(m^o,nr) values for the circular well, with 
a resulting revival behavior consistent with Tr^'^^\ since no m = states are allowed. 

Another variant would be an annular circular billiard, with an inner infinite wall at 
r = Rinner < R- The energy eigenstates can also be derived using Bessel function solutions 
(now including the 'irregular' or divergent Y\m\{kr) terms since the particle is kept explicitly 
away from the origin by the inner wall) with the energy eigenvalues resulting from the 
condition 

WKB-type expansions for the quantized energies are also useful in this case. The qualita- 
tively new features present in this geometry include new classical orbits (bouncing off the 
inner wall) but also diffraction features, as seen in periodic orbit theory analyses |2^] of such 
systems. (For the WKB analysis of the energies corresponding to orbits which bounce off 
the 'hard' inner wall, one must use (n^ + l) in place of (nj. + 3/4) in Eqn. While m = 

states are allowed, the special central, zero- momentum initial state is not, and whether a 
pattern of something like the Tr'^'^'^^ revivals is supported is currently under study. 
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APPENDIX A 

We briefly consider the special case of a Gaussian wave packet (of the form in Eqn. ()33p ) 
with vanishing momentum and initially located at the center of a two-dimensional square 
infinite well (or billiard) of dimension L x L. Because the problem is entirely separable, the 
autocorrelation function for the 2D problem will be a product of the individual ID values, 
namely A{t) = Aj.{t) x Ay{t), so it suffices to consider the ID case. The energy eigenstates 
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and eigenvalues are given 

Un{x) = —sin {^-j-j and = ^^^^ with n = 1, 2, 3, 4, ... (Al) 

(over the range (0,L)) and the general revival time is given by Trev = ^mL'^/hTi. For an 
initial Gaussian with vanishing momentum (pox = 0) and located at the center of the well 
(xo = the ID expansion coefficients, a„, simplify since the 'odd' parity states (here 

meaning n = 2,4,...) make no contributions to the wave packet and the energies can be 



written in the form 



E^ = ^-^{2n-lY with n = 1,2,3,... (A2) 



or as 



= 1^ [4n^ -^n+l]=^^ [8n{n - l)/2 + 1] = 8 [N{n) + l/s] (A3) 

where N{n) = n{n — l)/2 is an integer (neither even nor odd in general). Just as with the 
m = case of the circular well, in this very special alignment, the modulus of A{t) is unity 
with a reduced revival time of TrlT^^^^ = Trev/8, with a predictable phase factor (due to the 
constant, 1/8 term) at integral multiples of Trev^^^'^\ If one moves away from this special 
case by having xq L/2 or 7^ 0, this special revival structure is lost and only the (still 
exact) Tj-ev revivals are evident. 

We illustrate this in Fig. 4 for the case of a central, no-momentum solution (a) and for 
Xq ^ L/2 {h) and ~ (c) cases. We also show in the bottom (d) of Fig. 4 another special 
case where certain expansion coefficients vanish for symmetry reasons in the no-momentum 
case (xq = 2L/3, where every third a„ is zero) with exact revivals at integral multiples 
of Thus, a no-momentum 2D Gaussian wave packet moved slightly away from say 

(a;o,yo) = (2L/3,L/3) would experience the same kind of 'broken' revival time symmetry, 
as one moved from the center. 

A similar set of 'accidental' or 'symmetry' revival times exists for the equilateral triangle 
billiard. For example, for a zero-momentum state placed at the geometrical center, the 
revival time is Trev/Q where Trev is the exact revival time in Eqn. ((7j). Similar 'symmetry' 
points exist at distances of ^/3L/12 from the center in the direction of each vertex where 
the revival times are Not surprisingly, we find no additional such 'symmetry' points 

(besides the center) in the circular case. 
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FIG. 4: Plots of the one-dimensional autocorrelation function, A{t) versus t, for various Gaussian 
wave packets in a one-dimensional infinite well defined over the length (0,L). For case (a), where 
xq = L/2 and po = 0, there are special revivals at T^e^/S due to the extra symmetries forcing all 
even expansion coefficients to vanish. The values of |a„| are shown directly to the right, with arrows 
indicating those which vanish identically for symmetry reasons. Cases (b) and (c) respectively show 
the effect of changing xo and po slightly away from the values in case (a), illustrating how only 
the 'true' revival time is maintained. In case (d), we show another special case {xq = L/3) where 
certain a„ vanish (every third one in this case) for symmetry reasons, also giving accidental revival 
times. 

APPENDIX B 



While we have focused on the longer-term, revival dependence of wave packets in the 
circular well, it is interesting to note how the information about the classical closed (or 
periodic) orbits in this system is encoded in the energy eigenvalue spectrum, especially 
since most of the experimentally observed 2D circular billiard systems Q|, have made 
measurements which are relevant for short-term, quasi-classical ballis tic p ropagation. Such 
closed orbits are also the ones of relevance to periodic orbit theory 17[ measurements of 
such billiard systems. 

n 

For a system with two quantum numbers, there are two classical periods which in 
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our case are given by 



rpin,.) _ 

'''cl — 



2TTh 



and 



T. 



(m) 



(Bl) 



dE/dur ^' dE/dm 

and the two periods can beat against each other to produce the classical periodicity (T^°) 
for closed or periodic orbits if they satisfy 



T. 



(po) 



d 



(B2) 



with p > 2q for this geometry. We can then use this formalism to understand how these 
conditions can give rise to the classical expressions for the minimum radius and path lengths 
in Eqn. (j50|) . Instead of using the approximate expression in Eqn. (j46|) for the (m, -n-r)- 
dependent energies, we make use of the WKB condition in Eqn. (jlfij) and simply take partial 
derivatives of both sides with respect to n,. and m respectively. We thus obtain the conditions 



2h? 



dr 



2h^ 



R 



R^,„ \/E- m?'h?/2^r'^ 
dr ( dE 



dE 
drir 



'E -m?h^/2fir^ \dnr /ir^ 
The condition to be satisfied for periodic orbits can then be written as 



(B3) 
(B4) 



q _ T^i^^ _ {dE/dm) _ ( \m\h \ 
p ~ ~ (dE/drir) ~ \7Ty/2jlE) 



R 



dr 



R„ 



r^/r 



2 

iriin 



(B5) 



Evaluating the integral and using Rmin = \m\h/y/2jlE, we find that 

q 1 -if R 

- = — sec 

p 71 \Rrm 



or R^in{p,q)=Rn 



it COS I 

p 



(B6) 



as the condition on periodic orbits, as expected. To find the classical period for such closed 
orbits, we note that 



T, 



(po) 



27ihp 



2p\ R^ - R 



2 

"min 



/i [2pRsui{-Kq/p)] L{p,q) 



(B7) 



(dE/dnr) V~" ^ " "™"7 V 2^ 

where we identify Vq = ^/2EJJ1 with the classical speed. 

The classical periods for the closed orbits for the 2D annular well mentioned above can 
also be obtained from the WKB approximation in the same way, including those for the 
new features corresponding to 'bounces' off the inner infinite wall when Rmin is replaced by 
Rinner the integrations. The classical periods for the closed orbits for the 2D square and 
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equilateral triangle billiards can, of course, also be obtained in the identical manner, using 
the exact results for their energies. 
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